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Abstract 

We study the partition function of the model formulated with Wil- 
son fermions with only one species, both analytically and numerically. 
At strong coupling we construct the solution for lattice size up to 8 x 8, 
a polynomial in the hopping parameter up to 0(k^^^). At /? > we 
evaluate the expectation value of the fermion determinant for com- 
plex values of k. From the Lee- Yang zeroes we find support for the 
existence of a line of phase transitions from (/3 = 0, k ~ 0.38) up to 
(p = oo,K = 1/4). 
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1 Introduction 



QED2, the theory of electrons and photons in 2D, for massless electrons is 
analytically solvable ^ and has been studied extensively 0. In the original 
version the system has n/ = 1 fermions. The fermions are confined and the 
model is equivalent to a system of non-interacting bosons. The theory is 
superrenormalizable and there is only finite renormalization of the charge. 
Chiral symmetry is broken due to an anomaly in the axial current thus there is 
no Goldstone boson, but a massive pseudoscalar of Gaussian nature. One has 
charge shielding, i.e. there are no long-range forces between static charges. 
Allowing for Uf > 1 fermions, one does have additional massless states |Q, a 
situation prototyped in the 0{nf) non-linear a-model. 

Many of these properties are intriguingly similar to those of QCD in 
4D, which in the non-perturbative domain is studied mainly in the lattice 
regularization by computer simulations. The Schwinger model formulated 
on a lattice is therefore a challenging 2D model for lattice QCD4. In this 
formulation it naturally encompasses the massless and the massive situation. 
In the continuum the massive Schwinger model cannot be solved explicitly 
although there are results perturbative in m/e |Q, |^, |^. From these we expect 
a situation similar to the massless case with further bound states and quark 
trapping, again very much like in QCD4. 

Lattice formulations of fermions are plagued with the doubling problem. 
There are more than one pole in the Brillouin zone of the momentum space 
propagator. In the staggered (Kogut-Susskind) formulation one distributes 
different components of the spinor on different sites of a hypercube and 
thereby effectively reduces the 2^ multiplicity by a factor 2[^/^l. In the 
Wilson formulation the doubler modes are given masses of 0(1/ a) and one 
expects decoupling in the continuum limit a 0. Furthermore most lat- 
tice calculations rely on the hybrid Monte Carlo method to implement the 
fermions. In this process, however, another doubling of fermion modes is 
introduced in order to deal with the positive definite square of the determi- 
nant instead of the determinant itself. The only known way to avoid this 
doubling is to include the determinant in the observables and not use it as a 
probability weight. Thus, in the standard lattice formulation one deals with 
an Uf > 1 Schwinger model. Recently some exact results became available 
for the continuum massless nj > 1 Schwinger model 0, thus the situation 
in respect to the flavour problem of lattice versions has somewhat improved. 
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The naive and the staggered formulation has a simple phase diagram in 
the m)-plane, where (3 = 1/e^ is the gauge coupling and m denotes the 
bare fermion mass parameter. The continuum limit is recovered approaching 
the point (oo, 0) along trajectories with m ~ where the proportionahty 

constant characterizes the physical mass of the bound state boson. In that 
limit the known condensate value -{ipip) should be obtained. Most of the 
lattice calculations worked in that formulation and reproduced the expected 
continuum results of the massless Schwinger model successfully (see e.g. ^ 
1^, for recent work on the non-compact formulation cf. 0]). 

Little is known for the Wilson formulation, though. Here the two bare 
coupling parameters are the gauge coupling and the hopping parameter k 
related to the bare fermion mass. Chiral symmetry is explicitly violated 
for finite f3 and most likely restored in spontaneously broken form in the 
continuum limit. This deprives us of a suitable order parameter. Also, for 
the Uf = 1 model, there is no massless Goldstone boson. A continuum limit 
should be obtained in the approach to (/3 = cxo, k = 1/4) along trajectories 
following fixed values of physical (dimensionless) quantities. In QCD^ there 
is a line ndP) where the pion mass vanishes, corresponding to the situation 
of a vanishing bare fermion mass (due to the Wilson term the bare quark 
mass is renormalized additively). In the Schwinger model we may expect a 
similar behaviour and in this case the massless continuum Schwinger model 
is obtained following such a curve n = ndP)- 

Despite good knowledge of the properties of the theory in the continuum, 
it is still a challenge to clarify the nj = 1 model phase diagram on the lattice 
in the Wilson formulation. One cannot use (tp-ip) as an order parameter nor 
can we determine the mass of the boson state in a direct simulation. In this 
paper we therefore try to contribute to this issue in an analytic (at /9 = 
and oo) study and a computer calculation (at /3 = 1 and 5) of the partition 
function itself. With help of the equivalence of the strong coupling model 
to the 7-vertex model PU| , |TT[| we construct the hopping expansion up to 
0(k^^^). At non-zero (3 we directly determine the fermion determinant in 
the numeric integration over the gauge field background. The scaling of the 



Lee- Yang zeroes indicates that, somewhat contrary to an earlier finding [11 



there is indeed a phase transition line for all positive values of the gauge 
coupling /3. A subset of the numeric results has been presented elsewhere 



3 



2 Action, partition function and its zeroes 



The action for the massive lattice Schwinger model in the Wilson represen- 
tation is given by 

S{^,P) = Sf{k)+PSg, (1) 

+V'(a;)(l - -iy)U^{x)i){x + jl)) - ij{x)ilj{x)^ 
= J2 ^{x)[Ma{k,U)U^>^{x') . (2) 

For the gauge fields So is the standard Wilson action with the lattice gauge 
coupling P and Ma(k, U) denotes the Dirac operator on the lattice A. Eq.s 
(|l]) and are formal expressions in the partition function of the model, 
which, after Grassmann 'integration', is given by 

ZA{K,f3) = y"d/i(f/)detMA(/t,t/)e-^^G(^) (3) 
= Zc,,A(/3)(detMA(«:,f/))G. (4) 

^g,a(/3) is the partition function of the purely bosonic system and {0{U))c 
the expectation value of some operator for that system. Without loss of 
generality we normalize Zg,a{^) = 1- The partition function Z\{k) is a 
polynomial of degree 2|A| in /t and thus an entire function. All coefficients 
of that polynomial are positive. We know that det Ma is strictly positive for 
n{[3) < Kc{(3 = oo) = 1/4; above this value of k we checked explicitly that 
individual configurations do produce negative values although the partition 
function remains positive for finite lattices. 

The zeroes of that polynomial, called Lee- Yang zeroes [0, have a non 



vanishing imaginary part for all finite |A|. In the thermodynamic limit these 
zeroes pinch the real axis of the complex k plane and define thereby the point 
of non-analyticity at Hc- How the zeroes approach the real axis is governed 
by finite size scaling, at a critical point related to the critical exponents of 



the system [14 



From we see that the partition function of the full model is pro- 
portional to the purely bosonic expectation value of the determinant of the 
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lattice Dirac operator. Since Zg,k{I3) > we find that Z\{n^l3) is zero only 
where (det Ma(k, t/))G = 0. Whereas mass gap calculations cannot be per- 
formed for one flavour, this operator can be calculated and this calculation 
can be extended to any number of flavours. 

For given /3 we have to determine the expectation value for the deter- 
minant for complex k G C. Only at = and (3 = oo this may be done 
analytically. Elsewhere one may rely on the following numerical methods. 

analytic continuation: One obtains values of (det Ma) for various real k. 
Performing a polynomial fit and using this fit one analytically continues 
to complex k. In a test for the controllable situation of free fermions 
we find unstable results; even the closest zeroes could not be identified 
reliably. A better alternative is to obtain the coefficients of the poly- 
nomial expressed in terms of moments by direct numerical simulation. 



This approach has been successfully used in a recent study in 4D [|15 
The crucial point is the convergence of the coefficients with n. This 
approach is closest related to the quite successful method of analytic 
extrapolation and determination of zeroes via histograms [jl6| . 



direct evaluation: One integrates over the gauge fields with the usual 
Monte Carlo simulation and produces gauge configurations with the 
standard bosonic measure e~^^°^^'>djji{U). For each configuration one 
determines det Ma for a sufficiently dense set of complex k in the ex- 
pected region of the closest zero. The gauge field average eventually 
produces the required results. This method is applicable for small lat- 
tices only, since one has to determine the determinant for each gauge 
field configuration explicitly. On the other hand it does not involve any 
extrapolation and works for arbitrary number of flavours. 

3 Analytic results 

3.1 Strong coupling limit (/? = 0) 

For small lattice size (e.g. 2x2) one may analytically solve 

2|A| 

detMA(«:,f/) = Y,Cn{UK (5) 

n=0 
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and integrate over the gauge fields explicitly to obtain 



Za(/€) = (detMA(/€,?7))G (6) 

as a polynomial in n. Although done with symbolic programs this becomes 
prohibitive for larger lattices. 

Another method is based on a map of the massive lattice Schwinger model 
(at /3 = 0) on the eight- vertex model |T^. Integrating over the gauge fields 



first one finds a theory of non-intersecting loops represented by a vertex 
model, 

ZsvAM) = ^ n ' (7) 

where at denotes the coupling for vertex type i and rii the multiplicity. The 
Schwinger model corresponds to ai = M^, 0^2 = 0, 0:3 = 04 = 1, aj>4 = ^ and 
M = 1/{2k). In particular the vertextype 2 corresponding to the intersection 
of loops is forbidden. We then produce all possible configurations summing 
up the corresponding weights. 

A configuration may be represented by a legal set of vertices or, equiv- 
alently, by a collection of non-intersecting closed loops (out of connected 
links). For the determination of the series we simultaneously use the link- 
and the vertex representation. For lattices of size L x N (periodic b.c.) the 
number of a priori possible link configurations is 2^^^. 

The transfer matrix approach is by far the most economic approach. For a 
lattice of size LxN we consider the transfer matrix for a column of L vertices, 
T(a, b) where a, b denote the link configurations of the left and righthand set 
of L links each, e.g. a = (ai, 02, . . . , a^); a and b may assume 2^ values each. 
The internal variables have been integrated and periodic b.c. conditions in 
the vertical direction have been implemented. Each entry in T is given by a 
polynomial term M^"^ if it has ni vertices of type 1. 

One has the symmetries T{a,b) = T{b,a) = T{d,b) = T{b,d), with the 
notation d = (a^, . . . , 02, cti), reversal of the direction. Furthermore only en- 
tries with even number of occupied external links (corresponding to values 
of a and b) and allowed vertex configurations are non-vanishing. This effec- 
tively reduces the number of relevant entries by a factor of roughly 8. In the 
computer implementation only relevant entries of T are kept. For L = 8, T 
has only 8384 independent non-zero entries, each a polynomial in M. From 
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T one constructs 

T\a,c)=J2T{a,b)T{b,c) , = T\T\ (8) 

6 

and so on. Then 

Zsv{M; LxN) = Tt{T^} = ^ T^(a, a) . (9) 

a 

Alternatively, we also considered blocks of 4 x 4, which may be described 
by functions, P{a,b,c,d), with four arguments denoting the link configura- 
tions of the four edges. Again one has various rotational and reflection sym- 
metries reducing the number of independent non-zero entries of P, which are 
constructed explicitly at the begin. From P(a, b, c, d) one may construct the 
series for e.g. 4x8 and 8x8 lattices by suitable summations, like 

Zsv{M;8 X 8) = ^ P{a,b,c,d)P{a,u,c,t)P{r,b, s,d)P{r,u, s,t) . (10) 

abcdrstu 

This method is less efficient than the one-column transfer matrix approach, 
which we finally used to construct the series. 

In the 7-vertex model the partition function is given as a series in the 
variable M, Zsv{M;L x N) = Y.ndnM''. In the 2D Schwinger model one 
uses the hopping parameter k and in table |I| we give the coefficients of the 
series in this variable, 

Z^{k) = (2«:)^^Z8y(T^; LxN)=Y. c„«" , (11) 

for various lattice sizes. The coefficients integers and are related to 

the coefficient of the series in M through c„ = T^diN-n- Fig. la give the 
positions of the zeroes in the first quadrant of C for lattice size 8x8. 

3.2 Free fermions (/3 = oo) 

For U{1) gauge systems with torus geometry the choice between periodic or 
antiperiodic boundary conditions (b.c.) for the lattice Dirac operator Ma 
is irrelevant for all beta. Assume a field configuration in fixed gauge, then 
antiperiodic b.c. amount only to multiplying the gauge fields at the boundary 
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with a factor of —1 G U{1). This new configuration is within the sum over 
all gauge field configurations with the same weight as the original one. Thus 
it is a symmetry of the gauge field integral. 

At j3 = oo, if we fix the gauge, the choice of b.c. may be parametrized 
by two U{1) group elements {A,B), e.g. for antiperiodic b.c. A = B = — 1. 
Thus, as is done usually in analytic calculations, we may discuss the system 
of free fermions with antiperiodic b.c. as one particular configuration of the 
(5 — OG limit. In this case the fermionic action can be diagonalized and the 
determinant evaluated by Fourier transformation. Within the limitations of 
available workstations this may be done with symbolic computer programs 
up to lattice sizes 32 x 32 and larger. However, the results for the complex 
positions of zeroes are quite different for periodic and antiperiodic b.c. al- 
though the scahng behaviour agrees. Due to the dependence on the b.c. we 
cannot compare the zeroes directly with our results at finite j3. 

Unlike < oo in the free case the zeroes occur degenerate with multiplic- 
ity 4 and 8; fig. lb shows the zeroes (in the first quadrant of C) at = cxd 
for lattice size 8x8. For both, at /3 = and /3 = cxd, we find a distribution 
of zeroes, that does not follow a simple geometric shape. The general shape 
of the distribution aX (5 — oo becomes clearer at larger lattices, as shown in 
fig. Ic for lattice size 32 x 32. 

4 Numerical results and discussion 

For < /3 < oo we use numeric techniques to obtain information on the 
closest zero on lattice of size 2 x 2, 4 x 4, and 8x8, following the second 
method discussed at the end of sec. 2. We simulated background gauge field 
configurations for f3 — 5, f3 = 1, and also for /3 = 0, in order to check 
the reliability of the approach. For each configuration we evaluated the 
determinant on a grid of 20 x 20 complex k values in the presumed region of 
the closest zeroes. The final expectation values on this grid are then analyzed 
with help of interpolation. For lattice 2x2 size we summed over 10000 gauge 
field configurations, for 4x4 and 8x8 over 5000. 

The analytic results at /5 = are in excellent agreement with the numeric 
results for lattice sizes 2x2 and 4x4. At the largest lattice size, however, the 
exact position is off from the numerically determined position by 2 standard 
deviations (the jackknife statistical error). Actually, P — is the worst 
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case in the gauge field integration, since the configurations are completely 
random. We have to conclude that for the largest lattice size the number 
of gauge configurations considered is too small. (A test at /? = with a 
statistics of 32000 did not sufficiently improve the situation.) For this reason 
for lattice size 8 x 8 we doubled the statistical errors obtained at /3 = to 
define the errorbars for the other /9-values and consider these values with 
some caution. 

Fig. 2 summarizes our results for the size- and /^-dependence of the zeroes 
closest to the real /t-axis. The imaginary parts show the tendency to vanish 
for L — s> oo indicating the existence of a phase transition for all f3 along 
a curve Kc{(3). For free fermions, where k — Kc — 1/^, one finds 0{1/L) 
dependence. For /5 < oo, for the lattices considered, the behaviour indicates 
an even faster approach towards zero, like the 0(1/L^) expected at first 
order transitions (where the susceptibility ~ L^). However, the errorbars and 
smallness of the lattices do not justify stronger statements. Also, considering 
the non-uniform behaviour of the real parts shows that we are not yet in the 
asymptotic regime. 

As mentioned, we know of no suitable order parameter to identify exis- 
tence and type of the possible phase transition line. Earlier investigations of 
the behaviour of the strongly coupled massive lattice Schwinger model based 
on the 7-vertex model []11| rather suggested that there is no second critical 
point at /5 = for K < oo. To clarify this inconsistency we repeated that 7- 
vertex calculation with significantly increased statistics (up to a factor of ten) 
and denser grid in the coupling. It turned out that with this improvement 
one does indeed find indications of scaling in the susceptibility for lattices 
larger then |A| = 32 x 32. There is still a controversy concerning the peak 
position and the boundary conditions of the 7-vertex model to be settled. 

In summary, the n/ = 1 massive (and massless) Schwinger model on the 
lattice and with Wilson fermions is still a formidable task. Direct simulation 
of the model and determination of the bosonic mass seem to be impossible 
with present day resources. From the equivalence to the 8-vertex model we 
know that there is an isolated critical point at /? = and k, = oo |TD|, [T^]; 
this is most likely not the endpoint of the singular line discussed above. 
The expansion in k, around (/5 = 0, k = 0) is convergent only for |k| < 
R < 1/2[|11[]. The responsible singularity could lie at any complex k, with 
\k\ = k and therefore k does not necessarily define a critical point. The 
model is believed to have a line of phase transitions at Kc{i3) beginning at 
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Kc(cxd) = l/[2d) which is a point with a second order transition and running 
to some fi;c(0) < oo. Except for /3 — >• oo, the free model, there is no proof 
that this statement is correct. 

Our findings from the direct evaluation of the partition function and the 
finite size dependence of the Lee- Yang zeroes support the scenario of a phase 
transition line from /? = 0, k ~ 0.38(2) up to = oo.k = 1/4. The phase 
diagram may look similar to the Uf > 1 situation, also for QCD4. 
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Table 1: Series for Za(«;; L x N) ^ EnCnK^- 



n 


2x2 


4x4 


8x8 





1 


1 


1 


4 


64 






8 


768 


2304 


1024 


12 




32768 


32768 


16 




2617344 


2547712 


20 




51904512 


121634816 


24 




1068498944 


6214909952 


28 




9663676416 


291286024192 


32 




37597741056 


16264921612288 


36 






733659871051776 


40 






32983549241982976 


AA 






1393668629299462144 


48 






57045024334275411968 


52 






2178933957104869834752 


56 






78035756256081760747520 


60 






2591089287867411030081536 


64 






79973512017556234532028416 


68 






2269245313085843057330356224 


72 






58685404081509064264407580672 


76 






1378038126803683051771720105984 


on 






zy i iUoo / 4zozzoooOoyUo4it)UzU4UUo4 


84 
88 






547924220534282839622026012917760 
9113040973297576534851687523811328 


92 






132393136404532724219119019631837184 


96 






1658720256337887142482109949475291136 


100 






17660710533168127842880532530400854016 


104 






157058824477979862076814317113632096256 


108 






1140500578260786292933495730676818771968 


112 






6579687246247908061970155944177873977344 


116 






28956384933789738373830226402627572203520 


120 






92150412611536628441916337226044921610240 


124 






187529160722546587292381268726888701886464 


128 






223167089080216837357579582404487997816832 
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Figures 



Fig 1: Distribution of the complex zeroes for (a) /3 = on an 8 x 8 lattice 
and at /3 = oo on (b) 8 X 8 and (c) 32 x 32 lattices. 

Fig 2: (a) The zeroes closest to the real K-axis at /9 = (squares), l3 = 1 
(crosses) and /3 = 5 (triangles); subsequent lattice sizes (2 x 2, 4 x 4, 8 x 8) are 
connected by lines to guide the eye. (b) Scaling behaviour of the imaginary 
parts of the closest Lee- Yang zeroes (notation as in (a)) vs. 1/L; we also 
show results for /3 = oo and antiperiodic b.c. as discussed in the text. 
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